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SECTION  I 
INTRODUCTION 

Gray,  Kelley,  and  Mclntire  CI9783  have  described  a  method  for 
determining  the  order  of  an  ARMA  process  and  for  identifying  roots  of 
the  characteristic  equation  on  or  near  the  unit  circle.  In  this  paper, 
we  will  demonstrate  how  that  approach  can  be  utilized  in  modeling  seasonal 
data . 

Many  practitioners,  at  present,  perfunctorily  employ  the 

S 

operator  1-B  on  any  data  set  believed  to  have  a  period  of  length  S. 

Use  of  that  operator,  however,  tacitly  assumes  not  only  a  frequency  of 

1/S  to  be  present  in  the  data,  but  of  all  the  harmonics  of  1/S  (i.e., 

it  assumes  the  frequencies  0,  1/S,  2/S,...,  Cs/2]  where  t.l  is  the 

S 

greatest  integer  function).  We  shall  demonstrate  a  technique  which  will 
aid  in  determining  when  such  an  operator  is  called  for  and  when  other 
seasonal  models  are  called  for. 

In  the  next  section,  we  define  our  terms  and  give  the  theorems 
which  are  necessary  for  describing  the  method  which  we  employ.  In  the 
following  section,  we  illustrate  the  procedure  using  two  example  series: 
the  International  Airline  series  given  in  Box  and  Jenkins  □  976]  and  the 
so-called  Radio  series  given  by  Siddiqui  D962J.  >  \  I  I  T 

I  I  1  < 
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SECTION  II 

DEFINITIONS  AND  THEOREMS 
Definition  1 


By  an  ARMA  (p,q)  process,  we  mean  a  stochastic  process  {Xt>  which 
satisfies 

0(B)  Xt  =  0(B)at,  t  =  0,  +1,  +2,... 

where 

0(B)  «  1  -  0^  -  02B2-...-  0pBP, 

0(B)  =  1  -  0^  -  02B2-...-  0qBq  with 
W°  andBkXt  =  Xt_k 

The  algebraic  equation  0(r)  =  0  is  called  the  characteristic  equation 
of  the  corresponding  ARMA  process.  We  assume  that  0(r)  has  all  of  Its 
roots  on  or  outside  the  unit  circle,  and  that  0(r)  and  0(r)  are 
relatively  prime.  fat>  is  assumed  to  be  a  white  noise  process. 

It  is  well  known  that  {Xt>  is  a  stationary  process  If,  and 
only  if  all  of  the  roots  of  its  characteristic  equation  are  strictly 
outside  the  unit  circle.  By  a  non-stationary  ARMA  process,  we  will  mean 
an  ARMA  process  with  one  or  more  of  the  roots  of  0(r)  lying  on  the  unit 
circle.  That  is,  we  exclude  the  case  of  roots  Inside  the  unit  circle. 
The  term  seasonal  model  will  be  used  to  designate  the  following  class 
of  non-stationary  ARMA  processes. 

Definition  2 

A  factor  y(B)  will  be  said  to  be  seasonal  if  ( 

f(B)  *  1  +  B,  or 
?(B)  *  1  +  B  +  B2,  |?j|  <  2. 
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An  ARMA  (p,q)  process  will  then  be  referred  to  as  seasonal  if  it 
consists  of  one  or  more  seasonal  factors. 

Motivation  for  the  above  definition  Is  seen  most  easily  in  the 
frequency  domain.  A  seasonal  factor  is  any  (irreducible)  non- stationary 
factor  with  associated  frequency  greater  than  0. 

Definition  3 

The  autocorrelation  of  a  stationary  ARMA  process  is  given  by 
p(k)  =  E(XtXt+k)  /  E(X*). 

Strictly  speaking  the  autocorrelation  of  a  non-stationary  ARMA 
process  does  not  exist.  However,  if  one  regards  a  non-stationary  ARMA 
process  as  a  limiting  case  of  a  sequence  of  stationary  processes,  the 
following  definition  will  appear  natural. 

Definition  4 

Let  pk  (X-|,...,Xp,  £)  denote  the  autocorrelation  at  lag  k  of 

a  stationary  ARMA  process  with  ,  the  roots  of  <j)(r)  and 

9>s*  (8-j .... ,8q)  the  moving  average  parameters.  Now  suppose  that 

{Xj.}  is  a  non-stationary  ARMA  (p,q)  process  with  roots  of  0(r) 

Xt . x_  of  which  X, ,...,xm  are  on  the  unit  circle.  We  then  extend 

l  p  l  m 

the  definition  of  p(k)  by  letting 

p(k)  =  1  ini  pj^(ttA'i  fCiXg»  ♦  # »  » » *  ♦  >*p>8) 

a-l+ 

We  will  call  p ( k)  the  autocorrelation  of  the  non-stationary  process 
{Xt}. 
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Definition  5 

If  X-j,...,Xn  are  consecutive  random  variables  from  the  ARMA 
process  {Xt>,  we  take  as  our  estimator  of  p(k) 


Fundamental  to  the  discussion  to  follow  is  the  so-called  S-array. 
For  completeness,  we  give  a  formal  definition  of  the  S-array;  however,  its 
importance  for  our  purposes  is  contained  in  the  two  theorems  which  follow. 
(Our  definition  differs  slightly  from  that  normally  given.  It  represents 
a  simple  shift  in  Index  in  order  to  give  the  format  of  the  shifted  S -array 
suggested  by  Woodward  and  Gray  [1979]  in  simpler  notation). 


Definition  6 

Given  a  doubly  Infinite  sequence  {fm>  let 


1 

1 

...  1 

^m-n+1 

1 

^m-n+2 

1 

•••  firH-l 

1 

/ 

^m-n+1 

1 

•  •  • 

i 

1 

• 

1 

/ 

1 

i 

fm 

fmfl 

“•  fmfn 

/ 

Itl 

fm+n-l 

The  S-array  is  then  the  numbers  Sn(fm)  displayed  as  in  Table  1. 


Theorem  1 
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If  {Xf}  is  a  stationary  ARMA  (p,q)  process  and  if  f  =  p(m)  or 
t  m 

fm  *  (-1  )mp(m) ,  then  Sn(fm)  =  C-j  for  all  and  Sn(fm)  =  Cg  for  all 

mi.m2  and  only  if  n  =  p,  m^  =  q,  and  =  -q-1,  where  C.|  and  Cg  are 
constants. 

PROOF  SEE  GKM  [1978] 

We  thus  note  that,  given  the  true  autocorrelation  function  of  a 
stationary  ARMA  process,  the  S-array  provides  an  unequivocal  identifica¬ 
tion  of  p  and  q.  Given  an  estimate  p(k),  we  then  look  for  a  similar 
pattern  in  the  S-array  to  provide  information  as  to  the  order  of  the 
process  {Xt>. 

Theorem  2 

If  (X . >  is  an  ARMA  (p,q)  process  and  if  f  =  p(m)  or  if  f  =  (-1 )mP (m) 

t  m  m 

then  (Xj.)  is  non-stationary  if,  and  only  if  for  some  n  and  some  C, 

Sn(fm)  s  C  f°r  all  ni,  where  C  is  a  constant. 

In  that  case  n  is  the  number  of  roots  of  highest  multiplicity  among  those 
roots  of  <J)(r)  located  on  the  unit  circle. 

PROOF  See  Quinn  [1980]  and  Theorem  1  GKM  [1978] 

Theorem  2  suggests  that  a  stepwise  procedure  will  be  required  for 
identifying  the  complete  model  whenever  (j)(r)  has  roots  on  the  unit 
circle.  First,  the  presence  of  non-stationarities  are  detected  by 
noting  a  column  of  the  S-array  which  is  relatively  constant.  The 
series  is  then  transformed  by  the  indicated  non-stationary  factor  and 
the  residual  series  is  then  investigated. 
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Few  distributional  properties  of  the  S-array  are  known.  However, 

Gray,  Kelley,  and  Mclntire  [1978]  have  shown  through  a  variety  of 
examples  that  the  S-array  is  relatively  robust  to  stochastic 
disturbance.  We  also  give  two  asympotic  results. 

Theorem  3 

Suppose  that  are  consecutive  random  variables  from  an  ARMA 

(p,q)  process  satisfying  $(B)Xt=9(B)at  and  r^ ( k)  is  the  sample  autocorrelation 
function  at  lag  k. 

(1)  Suppose  Sn(p(k))  is  defined  and  P  -  lim  rT(k)  *  p(k),  then 

piHm  VrT(k,)  =  s><k»  T~ 

if  the  roots  of  <J)(r)  are  strictly  outside  the  unit  circle. 


( i 1 )  Xt  is  non-stationary  if,  and  only  if  for  some  n  and  C 

P  -  11m  Sn(rT(k))  =  C 
T— 

for  all  k,  where  C  is  independent  of  k,  and  Sn(p(k) )  is  defined. 
PROOF 

(i)  This  easily  follows  since  Sn  is  a  continuous  function 

and  P  -  lim  ry(k)  =  p(k). 

T-*» 

(ii)  The  proof  of  part  (ii)  relies  on  two  quite  useful 
results  which  were  established  by  Findley  [1980]  and  which  we  state 
as  lemmas. 
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Lemma  1 

Suppose  all  quantities  are  as  defined  in  Theorem  3  and  that 

m  A 

t(r)  =  [  ir(l-a,  r)]%(r) 

1  1 

where  the  a.  are  distinct,  |c^|  =1  and  ?{r)  has  no  roots  on  the  unit 

circle  of  multiplicity  greater  than  d  -  1.  Writing 
m 

n  (1-a.r)  »  1  +  a-,r  +  ...  +  a 
1 '  l  I  mm 

we  then  have 

P  -  lim  ^rT(k)  +  a-j r( k-1 )  +  ...  +  amry(k“m)^a  0 
for  all  k  *  0,  +1 . .  . 


!:.§?!!!§ If  we  let: 


rT(k-m+l)  .  .  .  rT(k) 


rT(k) 


rT(k+m-l ) 


the  denominator  quantity  in  the  S-function,  we  have,  taking  m  as  in 
Lemma  1, 

11m  inf  |Hm(rT(k))|  >  0 
for  all  k,  almost  surely. 

We  may  now  prove  part  (ii)  of  Theorem  3  using  the  notation 


introduced  above. 
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(=?)  .  Assume  Xt  is  non-stationary  and  0(r)  is  as  defined  in  lemma  7. 
Fix  k  and  let 

Uy(  i )  =  Tj(  k  +  i)  +  a^r^k+i-1 )+...+  a^  (k+i-m ) 
for  i  =  l,...,m. 

By  lemma  2,  P  -  lim  UT(i)  =  0,  i  =  l,...m. 

T-h®  1 

Further  by  performing  simple  column  operations  in  the  numerator 
determinant  we  have,  letting  A . .  be  the  ijth  cofactor  of  the  matrix 
in  the  numerator  of  the  S-function, 

sn/rT^k^  =  0  +  +  ...  +  am)  (-l)m  +  UT(1)  A2,nrH-l /Hm^rT^k^ ^ 

+...+  UT(m)  a„h.i  >nH.] 

HnWk» 

Now,  since  A^j  is  bounded  it  follows  from  lemmas  1  and  2  that: 

P  -  lim  Sm(rT(k))  *  (-l)m(l  +3^  ...  a  ) 

T-*» 

(^r)  Suppose  then  that  there  is  an  n  so  that 

P  -  lim  Sn(rT(k))  =  C  for  all  k 
T-*» 

and  that  Sn(p(k) )  is  defined  for  all  k. 

If  Xt  is  stationary,  then  with  the  above  conditions: 

P-  lint  Sn(rT(k})  =  Sn(p(k)) 

T"** 

which  is  not  the  same  for  all  k.  Hence  Xt  is  non-stationary  and  the 
theorem  is  proved. 


k 


J, 
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SECTION  III 
ANALYSIS 


To  illustrate  the  usefulness  of  the  S-array  as  a  model 
identification  tool,  we  consider  two  real  data  examples.  Our  first 
example  is  the  International  Airline  Series  given  by  Box  and  .lenkins 
[1976].  This  example  we  will  briefly  examine  even  though  it  is 
analyzed  in  much  the  same  manner  by  Gray  and  Woodward  [1980]  and  by 
Hart  and  Gray  [1980].  Our  purpose  for  including  the  Airline  Series 
is  to  show  the  contrast  between  it  and  our  second  example  series: 
the  so-called  Radio  Series,  given  by  Siddiqui  [  1962]. 

The  International  Airline  Series  consists  of  the  natural 
logarithm  of  the  number  of  passengers  in  International  air  travel. 

The  data  are  monthly  values  from  January  1949  to  December  1960  . 

A  plot  of  the  data  (see  figure  1)  shows  that  it  appears  to 
have  a  linear  trend  and  a  quite  distinctive  yearly  periodicity  about 
that  trend.  The  S-array,  evaluated  using  fm  *  (-l)mp(m),  is  shown  in 
Table  2. 

An  ambiquity  regarding  the  identification  is  noted,  since  both 
the  1st  and  13th  columns  are  relatively  constant.  The  series  will 
clearly  not  be  well -modeled  as  a  1st  order  process;  however  the  near 
constant  1st  column  indicates  the  presence  of  a  near  1st  order  non- 
stationarity.  It  Is  usually  best  to  remove  that  factor  before  attempting 
further  identification.  That  factor  is  estimated  as  roughly  1-.95B, 
using  the  Yule-Walker  estimate. 

We  next  transformed  by  the  operator  given  above.  The  S-array, 


Figure  1 

Airline  Data 


TABLE  2 
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using  f  =  p(m),  for  the  transformed  series  is  given  in  Table  3. 

There  is  now  an  unambiguous  constancy  in  the  12th  column.  The  Yule- 

Walker  estimate  of  that  operator  is  given  in  Table  4. 

At  this  point  a  subjective  decision  must  be  made  as  to  whether 

or  not  a  seasonal  model  is  desired.  If  a  seasonal  model  is  desired,  it 

must  be  decided  which  factors  to  alter  to  the  unit  circle.  From  the 

factors  given  in  Table  4,  it  is  apparent  that  each  of  the  frequencies 

12 

associated  with  the  operator  1-B  is  present.  We  also  note,  that, 
with  the  possible  exception  of  the  factor  1+.92B,  all  of  the  roots 
are  near  the  unit  circle. 

Thus  a  reasonable  approximation  to  the  2nd  estimated  operator 
12 

is  given  by  1-B  .  Likewise  the  operator  1-.958,  initially  estimated, 

may  be  adjusted  to  the  non-stationary  operator  1-B.  Other  possibilities 
might  be  considered  (see  Gray  and  Woodward  [1980]),  but  the  work  done 
so  far  indicates  that  the  operator  (1-B) (1-B  )  is  not  unreasonable. 
Table  5  gives  the  S-array,  using  fm  *  (-l)mp(m),  for  the 

1  p 

original  series  transformed  by  (1-B) (1-B  ).  The  residual  series 

appears  to  be  well-modeled  by  an  AR(12).  The  Yule-Walker  fit  is 

given  in  Table  6.  We  thus  arrive  at  the  model 

(l-B)(l-B12)0(B)Xt  =  at  0) 

where  {at>  is  a  white  noise  series  with  var  at  =  .00136  and  (>(B)  is  the 

stationary  operator  given  in  Table  6. 

12 

Taking  1-B  as  the  only  non-stationary  component  of  the  model. 


Gray  and  Woodward  p980]  arrive  at  the  model 
(1-B12)  $(B)Xt  *  at 


(2) 


TABLE  3 


ESTIMATED  WHITE  NOISE  VARIANCE  .001922 
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(12) fit  to  the  Airline  Series  after  being  transformed  by  1-.95B. 
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TABLE  5 


sn 

SI  2 

SI  3 

SI  4 

-63.460 

3.090 

-1.861 

.188 

-  1.330 

3.535 

-3.286 

-.270 

-  2.188 

3.514 

-41.970 

4.285 

-  8.356 

3.125 

-  1.442 

6.143 

1.075 

2.827 

-  9.870 

5.802 

1.364 

3.012 

-  .323 

2.478 

1.584 

3.067 

-  2.659 

2.781 

-  3.155 

3.069 

3.985 

10.676 

11.752 

3.282 

-11.955 

24.985 

-15.269 

2.812 

7.770 

11.889 

-  .711 

.952 

-  .848 

.914 

4.219 

1.152 

-  1.511 

1.152 

-5.292 

1.096 

-  .601 

-5.016 

.031 

1.064 

6.888 

-2.831 

-1.613 

1.057 

.125 

1.715 

-1.536 

1.036 

-1.491 

1.579 

-1.430 

1.064 

1.107 

1.825 

1.781 

1.292 

-1.396 

2.104 

-  .487 

1.280 

18.507 

-  .118 

-  .798 

1.401 

1.558 

.095 

27.383 

1.028 

.841 

2.240 

S-Array  for  the  Airline  Series  after  bging 
transformed  by  (1 -B) ( 1 -B ’ .  fm=(-l)mp(m) 


Estimated  White  Noise  Variance  .001364 
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where  0(B)  Is  a  13th  order  stationary  operator  and  {a^Jis  a  white  noise 
process  with  var  at  =  .001267. 

Gray  and  Woodward  [1979]  discuss  some  of  the  considerations  which 
are  relevant  in  deciding  between  models  (1)  and  (2).  They  argue  for  model 
(2)  based  on  the  models'  respective  forecast  functions.  They  also  give 
some  comparison  between  the  Box- Jenkins  model,  model  (2)  above,  and  a 
model  given  by  Parzen  [1979].  The  reader  is  directed  to  the  paper  for 
further  details. 

The  2nd  example  we  will  consider  is  the  so-called  Radio  Series 

given  by  Siddiqui  [1962].  The  data  consists  of  the  240  monthly  median 

^oF2  va^ues  observed  in  Washington,  D.C.  from  May  1934  until  April  1954. 

A  plot  of  the  series  and  of  the  autocorrelation  function  (see 

Figures  2  and  3,  respectively)  each  indicate  the  presence  of  a  low  frequency 

component  and  a  quite  distinctive  yearly  periodic  oscillation.  Noting  the 

12 

yearly  period,  many  practitioners  would  apply  the  transformation  1-B  to 
the  data.  Further  analysis  below,  however,  will  show  that  operator  to  be 
unnecessary  and  in  fact  deleterious  for  this  particular  data  set. 

The  S-array  using  fm=(-l )mi5(m)  is  given  in  Table  7.  The  1st 
column  is  seen  to  be  roughly  constant,  reflecting  the  low  frequency  component. 
Since  the  S-array  is  not  too  distinctive,  in  view  of  Theorem  2,  it  seems 
reasonable  at  this  point  to  prefilter  the  data  by  1-.9B.  At  this  point,  of 
course,  we  do  not  mean  to  imply  that  1-.9B  is  a  factor  in  the  model  but 
simply  an  appropriate  high  pass  filter  which  allows  clearer  identification 
of  the  model . 

The  S-array  for  the  transformed  series,  using  f  *  (S(m)  is  given  in 
Table  8.  The  S-array  is  relatively  constant  in  the  13th  column  which 

.  i 

indicates  that  the  untransformed  series  should  be  well  fit  by 


an  AR(14).  The  Yule-Walker  fit  is  shown  in  Table  9.  Note  the 


Figure  2.  Plot  of  the  first  236  points 
from  the  radio  series 


LAG  r(i«  MONTHS) 


Figure  3.  Plot  of  the  estimated  auto¬ 
correlation  function  of  the 
radio  series 
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TABLE  8 


sn 

S12 

S13 

SI  4 

-4.734 

.100 

1.118 

-2.119 

1.004 

-1.305 

1.181 

-2.482 

1.415 

-7.103 

.398 

-2.023 

-  .343 

.306 

1.539 

-2.217 

-  .296 

-17.485 

1.217 

-2.174 

2.186 

-1.718 

1.433 

-  .988 

-1.429 

-  .519 

1.906 

-1.586 

2.412 

-1.229 

1.866 

-59.006 

-  .988 

.548 

-.423 

.420 

.527 

.145 

-.407 

-2.300 

-  .543 

.296 

-.464 

-  .5  00 

.081 

-.092 

-.225 

.659 

.094 

-3.005 

-.290 

.643 

-  .364 

.348 

-.197 

.823 

-  .278 

-3.313 

-.661 

.787 

.764 

-  .069 

-.413 

.932 

-3.204 

.450 

-.441 

.422 

S-array  for  the  radio  series  after  being 
transformed  by  1-.9B. 


Estimated  White  Noise  Variance  .437271 
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AR(14)  Yule-Walker  fit  to  the  Radio  Series. 


factors  l-( .94+.05i )B  corresponding  to  the  low  frequency  component  of  the 
data.  Since  the  imaginary  part  is  small,  those  two  factors  are  close  to 
being  a  double  root  of  one.  Noting  Theorem  2  above,  it  is  not  surprising 
that,  with  only  an  estimate  of  p ( k)  available,  the  constant  behavior 
appears  in  the  first  column.  That  is  often  the  situation  when  there  is 
a  strong  low  frequency  oscillation  in  the  data.  It  should  be  noted, 
however,  that  an  oscillation  is  also  apparent  in  the  first  column  of  the 
S-array  which  is  contrary  to  the  typical  behavior  of  a  process  with  a 
double  root  of  one. 

Examination  of  the  factors  given  in  Table  9  indicates  that  not  all  of 

12 

the  frequencies  associated  with  the  operator  1-B  are  present.  It  also 

indicates  that  some  of  those  frequencies  which  are  present  are  too  far 

from  the  unit  circle  to  be  regarded  as  non-stationary.  Thus  the  operator 
12 

1-B  is  clearly  not  part  of  the  model  for  this  data  set.  In  fact  use  of 
that  operator  causes  the  low  frequency  component  to  be  mistaken  for  some¬ 
thing  very  near  to  a  double  root  of  one  (that  is,  the  model  which  results 

12 

by  first  operating  on  the  data  by  1-B  has  one  root  of  one  and  one  real 
root  slightly  larger  than  one).  That  causes  forecast  functions  for  the 
model  to  have  a  nearly  linear  component  to  them,  which  is  clearly  an 
unreasonable  forecast  function  for  this  data  set. 

As  with  the  last  example,  we  are  now  faced  with  the  subjective  decision 
as  to  whether  a  seasonal  model  is  desired  and,  if  so,  which  factors  should 
be  made  seasonal . 

The  factors  associated  with  the  frequencies  1/12,  1/6,  1/4, 
and  1/123*  seem  the  reasonable  candidates  to  be  made  into  seasonal 
factors.  However,  models  including  the  factor  associated  with  the 

*  Since  the  low  frequency  component  is  thought  to  be  associated  with 
sunspot  activity  we  have  given  it  the  period  estimated  in  Woodward 
and  Gray  [1978]  for  the  sunspot  series. 


24 


frequency  1/123  have  forecast  functions  which  are  very  unstable.  To 
show  the  effect  of  choosing  different  seasonal  models  we  consider 
the  two  models  obtained  taking  the  factors  associated  with  1/12  and 
with  1/12  and  1/6  as  seasonal. 

The  S-arrays  for  the  data  filtered  by  1-1.732B+B2  and  by 
2  2 

(1-1.732B+B  )(1-B+B  )  are  given  in  Tables  10  and  11,  respectively. 

2 

The  residual  after  operating  by  1-1.732B+B  appears  well-modeled  by 
an  AR(12),  as  would  be  expected  assuming  the  process  to  be  well- 
modeled  as  an  AR(14).  The  Yule-Walker  fit  is  shown  in  Table  12  and 
is  seen  to  be  very  similar  to  that  shown  in  Table  9  (disregarding 
the  factors  already  removed). 

Examination  of  Table  11  indicates  that  the  residual  from 
the  operator  ( 1 -1 . 732B+B2) ( 1 -B+B2)  is  well-modeled  as  an  AR(13).  The 
Yule-Walker  fit  to  that  model  is  given  in  Table  13.  The  higher  order 
indicated  here  may  be  a  consequence  of  Theorem  2  in  that  terms  which 
were  masked  before  removing  the  non-stationarities  are  now  apparent. 

Forecast  functions  of  various  lengths  were  calculated  from  a 
number  of  origins.  The  two  models  performed  similarly  in  the  cases 
considered.  Two  fairly  representative  forecast  functions  are  given  in 
Figures  4  and  5.  We  thus  have  two  quite  tenable  models,  both  of  which 
explain  the  data  very  well.  Which  model  to  use  will  depend  on  the 
uses  to  which  the  model  will  be  put,  and  on  what  if  any  physical 
significance  can  be  found  in  the  extra  parameters  which  were  fit. 

An  Important  point  of  the  above  example  is  that  the  operator 

$ 

1-B  should  not  be  used  indiscriminately.  The  radio  series  is  an 
example  for  which  a  cursory  examination  suggested  the  operator  1-B12  to  be 


TABLE  10 


S9 

S10 

Sll 

S12 

SI  3 

SI  4 

1.237 

-1.285 

.018 

.715 

-1.801 

3.812 

3.822 

.620 

-  .755 

.731 

-1.306 

.301 

.897 

.792 

-5.444 

.309 

-1.657 

1.113 

1.496 

.107 

-  .013 

1.135 

-1.624 

-4.289 

1.012 

.031 

-1.291 

1.150 

-1.901 

-1.943 

1.041 

.894 

-  .954 

1.056 

-  .895 

1.690 

.510 

.011 

-  .559 

1.180 

-3.655 

-3.696 

.872 

.541 

-  .554 

1.443 

7.162 

7.757 

-.120 

.154 

-  .213 

.250 

-  .241 

.249 

-.402 

.005 

-  .154 

.290 

-  .479 

-  .511 

-.237 

.151 

-  .153 

.324 

1.016 

.875 

-.327 

.010 

-  .031 

.183 

-  .412 

.282 

-.318 

.030 

.022 

.163 

-  .558 

.386 

-.455 

.158 

-  .161 

.168 

-  .533 

-1.160 

-.336 

.137 

2.793 

.434 

-  .568 

-  .126 

-.984 

-.498 

.014 

.250 

-  .421 

.640 

.019 

.991 

-  .256 

.254 

-  .338 

-  .270 

S-array  for  the  radio  series  after  being  transformed 
by  (1-1.732B+B2).  f^m). 


TABLE  11 


S10 

Sll 

SI  2 

SI  3 

SI  4 

.553 

.309 

2.400 

-.855 

.551 

-.015 

-3.177 

3.430 

-1.421 

1.585 

3.124 

-3.191 

-7.398 

-1.235 

4.456 

.565 

-1.079 

.814 

-1.595 

-2.654 

.596 

-1.272 

4.347 

-  .305 

.269 

.120 

-1.756 

-2.903 

-  .273 

16.598 

1.464 

-1.626 

-  .906 

-1.807 

1.892 

1.348 

-3.320 

.292 

-1.646 

-6.346 

1.178 

-1.163 

2.255 

-1.799 

1.118 

1.211 

-56.977 

-1.012 

-2.111 

4.165 

3.583 

3.490 

3.434 

-1.996 

-11.104 

.248 

-  .232 

.249 

-  .284 

.277 

.470 

-  .501 

-  .205 

-  .327 

.627 

.491 

11.729 

.483 

-  .309 

-  .348 

.151 

-  .234 

-  .078  • 

-  .212 

.165 

.139 

-  .814 

.175 

-  .208 

1.603 

.081 

-  .488 

.313 

-  .072 

.081 

.454 

-  .502 

.841 

-  .080 

-  4.589 

.266 

-  .595 

-1.060 

-  .615 

.551 

.284 

-  .345 

.301 

-  .883 

1.375 

.007 

-  .431 

6.170 

-  .994 

4.511 

.433 

-  .455 

-1.406 

-  .467 

-  .296 

S-array  for  the  radio  series  after  being  transformed 
by  (1-1.732B+BZ)(1-B+B2). 


Estimated  White  Noise  Variance  .443774 
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appropriate.  Further  examination,  however,  showed  that  some  of  the 

12 

frequencies  associated  with  1-B  are  not  present  in  the  data  and 

that  some  of  those  present  are  clearly  not  on  the  unit  circle.  In 

12 

fact  as  already  pointed  out  the  use  of  the  operator  1-B  on  the  radio 
series  causes  one  of  the  most  salient  points  of  the  data  set—the  low 
frequency  oscillation— to  be  mistaken  for  essentially  a  double  root 
of  one.  So  if,  for  instance,  the  purpose  of  the  analysis  was  to  in¬ 
vestigate  the  fitted  model  for  evidence  that  sunspot  activity  influences 
radio  transmission  that  evidence  has  been  badly  obscured  if  not  lost. 

ACKNOWLEDGEMENT :  I  wish  to  thank  Professor  Henry  Gray  for  the 

numerous  helpful  comments  he  has  made  on  earlier  drafts  of  this 
paper  and  for  the  many  discussions  we  have  had  on  this  and  related 
topics. 


32 


REFERENCES 


Box,  6.E.P.  and  Jenkins,  G.  M.  (1976),  Time  Series  Analysis: 
Forecasting  and  Control. 

Findley,  David,  (1980),  Personal  Communication. 

Gray,  H.  L.,  Kelley,  G.  D.,  Mclntire,  D.  (1978),  "A  New  Approach  to 
ARMA  Modeling",  Communications  in  Statistics,  B7,  pp.  1-77. 

Gray,  H.  L.,  and  Woodward,  W.  A.  (1980),  "Application  of  S-arrays  to 
Seasonal  Data",  Proceedings  of  the  1979  Tulsa  Time  Series 
Symposium. 

Hart,  J.  and  Gray,  H.  L.  (1980),  "Modeling  Seasonal  ARMA  Processes," 
ONR  Tech.  Report  No.  140,  SMU  Department  of  Statistics. 

Parzen,  E.  and  Pagano,  M.  (1970),  "An  Approach  to  Modeling  Seasonally 
Stationary  Time  Series,"  Journal  of  Econometrics. 

Quinn,  Berry  (1980),  "Limiting  Behavior  of  Autocorrelation  Function  of 
ARMA  Processes  as  Several  Roots  of  the  Characteristics  Equation 
Approach  to  the  Unit  Circle,"  Communications  in  Statistics,  B9, 
pp.  195-198. 

Siddiqui,  M.  M.  (1962),  "Some  Statistical  Theory  for  the  Analysis  of 
Radio  Propagation  Data,"  Journal  of  Research  NBS,  66D  (Radio 
Propagation),  pp.  571-580. 

Woodward,  W.  A.  and  Gray,  H.  L.  (1979),  "On  the  Relationship  Between 
the  R  and  S-Arrays  and  the  Box-Oenkins  Method  of  ARMA  Model 
Identification,"  ONR  Tech.  Report  No.  134,  SMU  Department  of 
Statistics. 

Woodward,  W.  A.  and  Gray,  H.  L.  (1978),  "New  ARMA  Models  for  Wolfer's 
Sunspot  Data,  "Communications  in  Statistics,  B7,  pp.  97-115. 


Unclassified 


SECURI TV  CLASSIFICATION  OF  THIS  FACE  (Whtn  Data  tntarad) 


REPORT  DOCUMENTATION  PAGE 


4.  TITLE  fantf  Subtitle) 

The  Determination  of  Seasonal  ARMA  Models 


7.  AUTHOR!  »; 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  RECIPIENT'S  CATALOG  NUMBER 


XL 


S.  TYRE  OF  REPORT  A  PERIOD  COVEREO 


Technical  Report 


C.  PERFORMING  ORG.  REPORT  NUMBER 

141 


(.  CONTRACT  OR  GRANT  NUMBER!*.) 


Michael  Morton 


N0001 4-75-C-0439 


S.  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Southern  Methodist  University 
Dallas,  Texas  75275 


II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Office  of  Naval  Research 
Arlington,  Va.  22217 


.  MONITORING  AGENCY  name  A  AOORSSSflf  HUtmnt  I 


I 


i  CmInIIIiii  Ottlca ) 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 


NR  042  280 


IZ.  REPORT  OATS 

October  1980 


I).  NUMBER  OF  PAGES 

32 


IS.  SECURITY  CLASS,  !•!  thlm  re port) 


la.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 


IS-  DISTRIBUTION  STATEMENT  fat  rAI*  Depart) 


This  document  has  been  approved  for  public  release  and  sale;  its 
distribution  is  unlimited.  Reproduction  in  whole  or  in  part  is 
permitted  for  any  purposes  of  the  United  States  Government. 


17.  DISTRIBUTION  STATEMENT  ft  tba  abatraci  antara*  In  Bleat  ZB,  It  IlfltrMl  tram  Depart) 


IS.  KEY  WOROS  (Centime  an  raearaa  cl*a  It  nacaaaatr  an*  Harttllp  I y  black  n umber) 


ZOMBOSTRACT  (Centtnue  an  raaaraa  al*a  It  naaaaaaar  an*  Identity  A f  black  mmbat) 


Gray,  Kelley,  and  Mclntire  (1978)  have  described  a  method  for  determining 
the  order  of  an  ARMA  process  and  for  identifying  roots  of  the  characteristic 
equation  on  or  near  the  unit  circle.  In  this  paper,  we  will  demonstrate 
how  that  approach  can  be  utilized  in  modeling  seasonal  data. 


oo ,sr„  M73  EOlTION  OF  '  NOV  SS  IS  OBSOLETE 
S/N  0 101*014*  AA0 1  l 


_ Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PAOE  fWiMi  Data  I 


